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Transition frequencies of Fe ll ion are known to be very sensitive to variation of the fine structure 
constant a. The resonance absorption lines of Fe ll from objects at cosmological distances are used 
in a search for the possible variation of a in cause of cosmic time. In this paper we calculated the 
dependence of the transition frequencies on (g-factors) for Fe ll ion. We found corrections to 
these coefficients from valence- valence and core-valence correlations and from the Breit interaction. 
Both the core-valence correlation and Breit corrections to the g-factors appeared to be larger than 
had been anticipated previously. Nevertheless our calculation confirms that the Fe ll absorption 
lines seen in quasar spectra have large g-factors of both signs and thus the ion Fe ll alone can be 
used in the search for the Q-variation at different cosmological epochs. 
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I. INTRODUCTION 

Recently discovered acceleration of the universe(for a 
review, see [l[) is usually regarded as evidence for the 
existence of dark energy. Cosmological evolution of the 
dark energy may cause variations of fundamental con- 
stants, such as the fine-structure constant a and proton 
to electron mass ratio fi — mp/nie (see, e.g., |2j] and ref- 
erences therein). Different models predict different be- 
havior of the coupling constants and it is extremely im- 
portant to measure a at different cosmological epochs^. 

Because of the relativistic effects in atoms their transi- 
tion frequencies depend on aZ, where Z is atomic num- 
ber. Therefore, one can study space-time variation of 
a by comparing atomic frequencies for distant objects 
in the universe with their laboratory values. To do this 
one needs to calculate so called g-factors for the atomic 
transitions which are defined as follows: 



wiab + qx, X = (a/aiab) 



1. 



(1) 



Generally g-factors may have different signs and scale 
with atomic number as Z'^ (see Q). That allows to look 
for a-variation by comparing transition frequencies of 
heavy and light ions 0, Q ■ However, under astrophys- 



ical conditions the relative fractions of different ions de- 
pend on the gas density fluctuations within the absorbing 
cloud and thus the measured positions of the absorption 
lines are affected by irregular Doppler velocity shifts of 
the bulk motions leading to an additional noise in the 
Aa/a measurements 0, 

Therefore, one can try to find an atom, where q-factors 
for different transitions have different signs and large ab- 
solute values. In Ref. [9| it was shown, that the line 
1608 A in Fe ii has large negative g-factor, while other UV 
resonance lines have large positive g-factors. Additional 
advantage of Fe ii is relatively small values of isotope 
shifts [S] . This fact was used in [1, [13, [HI, [13 to suggest 
the Single Ion Differential a Measurement (SID AM) pro- 
cedure which is less sensitive to uncertainties inherent to 
the method applied in ^ and allows us to measure 
l^a/a at a single redshift z with a sufficiently high accu- 
racy (^Ac/a 10"*^ - IQ-^). 

Here we report new calculations of the g-factors for 
Fe II with more accurate account for relativistic effects 
and electronic correlations. In particular, we studied cor- 
rections to the g-factors from the Breit interaction and 
the core-valence correlations. We found that both cor- 
rections are noticeable but do not change previous con- 
clusion 0] that Fe il has large g-factors of both signs. 



^ Each cosmological epoch is characterized by the redshift pa- 
rameter z, which is defined in the spectral observations as 
z = (Aobs - Aiab)/Aiab, with Aobs,Aiab being correspondingly 
the observational and laboratory wavelengths of an atomic tran- 
sition. Currently we observe objects with z ranging from (local 
universe) up to 7. 



II. METHOD OF CALCULATION 

In order to calculate g-factors one can solve atomic 
relativistic eigenvalue problem for different values of a. 
For example, if we calculate atomic frequency uj± for two 
values X = ±1/8 of the parameter x in Eq. ([T|), the cor- 
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responding g-factor is given by 

q^A{LO+-LO^). (2) 

The choice of a; = ±1/8 is usually sufhciently small to ne- 
glect nonlinear corrections and sufficiently large to make 
calculations numerically stable. The exceptions may oc- 
cur for the strongly interacting levels. Such cases require 
special consideration. 

Fe II ion has ground state configuration [Ar]3d^4s with 
7 electrons in open shells, hence it has rather dense and 
complicated spectrum, particularly in the astrophysically 
important region above 60000 cni^^ from the ground 
state. Calculations for Fe ii are very difficult because 
of the large number of valence electrons and high density 
of the spectrum [1, i, [H, [H . 

We use Dirac-Coulomb and Dirac-Coulomb-Brcit "no 
pair" Hamiltonians. At first stage we use the frozen-core 
approximation and solve valence eigenvalue problem by 
means of configuration interaction (CI) method. Using 
Eq. ([2]) it is important that CI space is the same for 
two values of x. When four-component finite basis sets 
are used to solve many-electron Dirac equations, the ba- 
sic orbitals depend on a. The incompleteness of the CI 
space can introduce some error for q-factors. To check if 
this error is small we make several CI calculations with 
different basis sets. 

At first we solve Dirac-Fock equations with the code 
to find core orbitals ls,...,3p3/2 and valence or- 
bitals 3(^3/2, 3(^5/2, 4s, 4pi/2 1 4p3/2- Then we add virtual 
orbitals, which are constructed using the method de- 
scribed in (itI . [isj . In this method an upper component 
of virtual orbitals is formed from the previous orbital of 
the same symmetry by multiplication by some monotone 
function of radial variable r. The lower component is 
then formed using kinetic balance condition. (An alter- 
native method to form basis sets was suggested in 19]). 
Our basis sets include s,p,d, and / orbitals with prin- 
ciple quantum number n < N and are designated as 
[Nspdf]. Configuration space is formed by single and 
double (SD) excitations from the configurations SdPAs, 
3d^Ap, and Sd^Asip. Respective dimensions of the CI 
space are significantly larger than in previous calcula- 
tions 

In order to study Breit correction to the q-factors we do 
calculations for Coulomb, Coulomb-Gaunt and Coulomb- 
Breit approximations, where Gaunt approximation ac- 
counts for the magnetic interaction and Breit interac- 
tion also includes retardation term. We have found that 
due to the large number of valence electrons and relative 
smallness of the nuclear charge Z it is necessary to in- 
clude Breit correction to the valence-valence interactions 
as well as to the core- valence and core-core ones. 

In our calculation the Ar-like core is rather small and 
rigid, with binding energies |e| > 3 a.u. For this rea- 
son one can expect that corrections from core excitations 
should be small. The same conclusion follows from the 
relatively good agreement between theoretical spectra in 
the frozen core approximation and the experiment [ssj . 



On the other hand, correlations between core electrons 
and valence 3d-electrons are enhanced. For this reason 
the core- valence correlations require thorough investiga- 
tion. 

Usually one can use the many-body perturbation the- 
ory (MBPT) to account for core-valence correlations. 
Unfortunately, it is difhcult to implement MBPT for 
atoms with open d-shells because valence d-shell electrons 
belong to the previous atomic shell. That leads to much 
stronger interaction between valence d-electrons and up- 
permost core electrons [13, Hi], [12]. Because of that it 
is impractical to use the V-^""" approximation, where 
MBPT is much simpler than in the case of V'^ potential 
with N ^ iVcoro- Generally MBPT includes a large class 
of the so-called subtraction diagrams, which cancel out 
for V^""'" approximation. That not only increases the 
complexity of calculations but also makes them less sta- 
ble because of the large cancelations between "normal" 
and subtraction diagrams [2S] . Due to these problems the 
CH-MBPT method [li|, which proved to be rather effec- 
tive for treating atoms with several s- and p-electrons in 
the open shells [H, H^, [13, [28, 29], is much less effective 
here. Because of that we combine CI-I-MBPT method 
with conventional CI calculation for 15 electrons includ- 
ing 8 outermost core electrons. 



III. RESULTS 
A. Valence CI 

The results of the 7-electron CI calculation of the 
spectrum of Fe ii are listed in Table HI The CI space 
for each basis set corresponds to SD excitations from 
configurations 3d^As, 3d^4p, and 3d^AsAp. It is conve- 
nient to divide all levels presented in Table U into three 
groups. The first group includes five levels of the even- 
parity ^Dj manifold. To the second group we relate the 
low-lying odd-parity states of the configuration 3d^Ap, 
namely ^-Dg^j 7/2' ^^11/29/2' ^^"^ ^^7/2- Finally, in the 
third group we include two levels of the configuration 
3d^4s4p (®-P°/2 ''^7/2) closely located states 

of the configuration 3d^4p: ^Gy2J ^-^7/2' ^"^^ '^-^7/2- 

Properties of the even-parity levels from the first group 
(fine structure splitting, g- and g-factors) are rather in- 
sensitive to configuration interaction. All these quan- 
tities change very weakly with increasing the basis set 
and the size of the configuration space. An explanation 
is that the weight of the leading configuration 3c?^4s for 
these states is about 95-97%. Small admixture of other 
configurations does not influence their properties. 

Behavior of the odd levels from the second group is 
rather similar. It is seen from Table [T] that g-factors of 
these levels remain practically the same while transition 
frequencies and g-factors change by less than 4% when 
the basis set is increased from [5spdf] to [Tspdf]. The 
weight of the leading configuration Sd^ip for the 
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TABLE L Seven-electron CI calculations of transition frequencies lo (in cm~^), g-factors, and g-factors (in cm^^) for different 
basis sets. Dirac-Coulomb Hamiltonian in the frozen-core approximation is used. Transition frequencies and g-factors are 
calculated in respect to the ground state. 
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"These g-factors do not include semi-empirical corrections from fitting g-factors. The respective corrections are discussed in 

Sec. Imp] 



state is about 90% and greater than 95% for other states. 
The levels of different multiplets with the same total an- 
gular momentum J are well separated and interact with 
each other rather weakly. 

Comparing calculated transition frequencies for this 
group of levels with experimental data we see very good 
agreement. Largest discrepancy does not exceed 3%. 
There is also a reasonable agreement between the q- 
factors obtained in this work and in Ref. @. Except for 
the state ScfAp^Py^ the differences between g-factors are 
at the level of 1-2%. For the state 3d^Ap ^P°i^ the differ- 
ence is about 5%. The explanation for this difference 
is probably the following. Calculation in the paper 
was done on the basis set [6spd/] and the CI space in- 
cluded SD excitations from the configuration 3c?^4p. In 
order to accelerate calculations this space was then trun- 
cated by excluding non-relativistic configurations, whose 
weight in the levels of interest was smaller than 0.5%. 
Here we do not truncate configuration space and also in- 
clude the SD excitations from the configuration 3d^4s4p, 
which corresponds to triple excitations from the configu- 
ration 3(i^4p. Admixture of the configuration 3(i^4s4p to 
the state 3(i^4p^P°y2 is at the level of 5%. For this rea- 
son the SD excitations from the configuration 3(i''^4s4p 
appear to be more important for this level. 

Table H] shows that transition frequencies for the third 
group of levels also rather weakly depend on the size of 
the CI space. Comparing the results obtained for the ba- 
sis sets [5spd/] and [Tspd/] we see that a most significant 
change of w (~ 320 cm^^) is for the ^Hy2 state. Even 
here it constitutes only 0.5% of the transition frequency. 
On the other hand, because of the high density of levels 



in this part of the spectrum, such changes are compara- 
ble to the spacings between different levels with the same 
J. The Lande-factors g for several levels with J = 7/2 
of the configuration id^Ap are changing rather strongly 
from one calculation to another. That signals that these 
levels interact with each other. However, this interaction 
does not affect the g-factors too much since all these lev- 
els have similar dependence on a. The Lande-factors of 
the levels of the configuration 2>d^AsAp, in contrast, are 
very stable. For the most interesting level ^Py2 theoret- 
ical (7-factor is somewhat larger than in the experiment. 
This will be discussed in detail in Sec. IIII Dl 

We can sum up, that on the stage of the valence CI 
calculation there is agreement on the level of few per- 
cent between theoretical and experimental frequencies. 
However, the theory does not reproduce the intervals be- 
tween levels in the higher part of the spectrum, where the 
density of states is very high. We see that theoretical fre- 
quencies do not change much with the size of the valence 
space. We can therefore conclude that we are rather close 
to saturation of the valence CI space and the remain- 
ing difference between theory and experiment is mostly 
caused by neglect of the core-valence correlations. We 
will discuss this part of the problem in Sec. IIII CI 



B. Breit corrections 

The results presented in Table U correspond to the 
Dirac-Coulomb approximation. As long as we are mostly 
interested in g-factors, which depend on the relativis- 
tic corrections to electronic Hamiltonian, the role of the 
Breit interaction is enhanced. In order to estimate corre- 
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spending corrections we made Dirac- Coulomb and Dirac- 
Coulomb-Breit calculations of the g-factors in the one- 
coniigurational approximation. The results of these cal- 
culations are given in Table [III 



TABLE IL One-configurational calculation of g-factors in 
Dirac- Coulomb, Dirac-Coulomb-Gaunt, and Dirac-Coulomb- 
Breit approximations (in cm~^); 5qgn and Sqbr are respective 
Gaunt and Breit corrections to Dirac-Coulomb g- factors qoc- 
The latter are calculated for x — ±1/8 using Eq. ([2j and also 
for X = ±0.01. 
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One can see that the Breit corrections to the g-factors 
are non negligible, particularly for the levels of the con- 
figuration 3d^4s4p. It is worth mentioning that, as usual, 
the magnetic part of the Breit interaction (Gaunt inter- 
action) is significantly larger than the retardation part. 
On the other hand, the Breit corrections to the valence- 
valence interactions are comparable to the corrections to 
core- valence interactions. The reason for that is rela- 
tively small size of the core and the large number of va- 
lence 3d-electrons. Because of that it is important to 
include Breit interaction not only in the mean field po- 
tential, but also in the residual two-electron interaction. 
On the other hand, the Breit corrections to the g-factors 
are sufficiently small (i.e. less than 10%), so that they 
can be calculated within one-configurational approxima- 
tion as in Table |TT1 Below we will add this correction to 
our final result in Table IIVI 

In Table |lT] we also checked the size of the nonlinear 
corrections from using Eq. ^ for evaluation of the q- 
factors. For the Dirac-Coulomb approximation we calcu- 
lated g-factors using x = ±1/8 and x — ±0.01. One can 
see that the differences between two calculations are on 
the order of a fraction of a percent. We conclude that 
nonlinear corrections to Eq. ^ can be neglected for all 
levels considered here. 



C. Core- valence correlations 

As we mentioned in Sec. Uthe size of the core- valence 
correlation corrections a priori is not clear. On the one 
hand the Ar-like core for Fe II is rather small and rigid. 
On the other hand valence 3d-electrons belong to the 
shell with the same principal quantum number as the 
outermost core electrons 3s, 3p. Thus, there is no spatial 
separation between core and valence electrons. 

There are two approaches that can be applied to ac- 
count for core- valence correlations. The first approach 
is to use the CI±MBPT method and another one 
is to include two 3 s- and six 3p-electrons into the va- 
lence space. The former method allows using long basis 
sets and accounting for correlations with all core shells. 
The main problem with this method is instability of the 
MBPT for the mean- field potential , which includes 
large number of valence electrons, i.e. for TV — TVcoro 3> 1- 
At present this method is implemented only for the sec- 
ond order MBPT corrections to valence Hamiltonian. 
That may be insufficient for accurate treatment of corre- 
lations of valence 3(i-electrons with 3s- and 3p-electrons 
from the core. It is known that second order MBPT usu- 
ally overestimates core-valence correlations. Potentially 
more accurate all-order methods, like the one suggested 
in Ref . [30| , are currently used only for the systems with 
two valence electrons. For example, coupled-cluster cal- 
culations of the g-factors of a number of divalent atoms 
and ions was recently reported in Ref. [Slj . 

The latter method is more straightforward, but the 
problem here is the enormous size of the CI space even 
for rather small basis set. Because of that it is impossible 
to saturate CI, which leads to significant underestimation 
of the core- valence correlations. Due to these problems 
neither method can be preferred and we used their com- 
bination to estimate corrections caused by core-valence 
correlations. 

We started with solving 7-electron CI problem for the 
basis set [5s4p(i]. Then we solved 15-electron CI prob- 
lem using the same basis set and found correction Sqcv,i, 
which accounted for the correlations between valence 
electrons and 3s- and 3p-electrons from the core. In both 
cases CI space included SD excitations. The dimension 
of the CI matrix for the 15-electron problem was on the 
order of lO''. That limited the length of the basis set 
we could use in this approach. Such a short basis set 
was obviously insufficient for accurate treatment of core- 
valence correlations. To account for the incompleteness 
of the basis set we performed two calculations of the q- 
factors using CI±MBPT method. In one calculation all 
summations over virtual states in MBPT diagrams ran 
over the same basis set as above and sums over core states 
included only 3s- and 3p-shells. 

Another calculation used a much longer basis set 
[21s 16p 21(i 19/ 14(7] and all core shells. Valence space 
for both calculations was the same as in the 7-electron 
CI above. The difference between these two calculations 
gave us the second correction Sqcv,2, which accounted for 
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correlations with the inner core shells and the incomplete- 
ness of the basis set used in 15-electron CI. The results 
of these calculations are presented in Table IIIII 

In the 15-electron CI calculation the odd-parity levels 
'^y2 ^^"^ ''^7/2 strongly interacted with each other and 
their g-factors significantly differed from the 7-electron 
calculation and from the experiment. That hampered 
comparison of the respective g-factors. We excluded 
these levels from Tabic Hill 

In general CI-I-MBPT even for a short basis set 
gave much larger corrections to g-factors and to tran- 
sition frequencies than 15-electron CI. Comparison of 
the transition frequencies with the experiment showed 
that CI+MBPT method strongly overestimated the core- 
valence correlations. On the other hand, the 15-electron 
CI was much closer to the experiment. Because of that 
we used CI-I-MBPT method only to account for the in- 
completeness of the basis set. Our final estimate of the 
correction from the core-valence correlation is given by 
the sum ^^cv = Sqcva + Sqcv,2- 



TABLE III: Core-valence corrections to g-factors, Sqcv ~ 
Sqcv,i + Sqcv,2 (in cm~^). The term 5qcv,i is a difference of 
15-electron and 7-electron CI on the basis set [5s4pd]. The 
second term accounts for the incompleteness of this basis set. 
It is obtained as a difference of 7-electron CI-I-MBPT on the 
basis sets [21s 16p 21d 19/ 14g] and [5s4pd]. 
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D. Level ^P7°/2 at 62172 cm" 



The level ^Py2 62172 cm ^ is the most important 
level for the search of a-variation. This one is the only 
level of the configuration 3(i^4s4p, which is observed in 
astrophysics. All other observed odd-parity levels belong 
to the configuration 3dP4p. Note that the ground state 
belongs to the configuration 3d^As. It means that corre- 
sponding transitions are of 4s — > 4p type, while transi- 
tion to the level ^Py2 is of 3d 4p type. The relativistic 
corrections increase binding energy of s-electrons and de- 
crease binding energy of d-electrons. Because of that the 



g-factor for ^Py2 ^Gvel is of the opposite sign to the q- 
factors of all other observed levels. The existence of the 
q-factors of both signs makes Fe ii ion so attractive for 
the search of a-variation. 

Initially there was some controversy around this level. 
In Moore's tables [s^l this level was erroneously marked 
as 3d^(JS)4p^Py2- It was an obvious misprint because 
there is no term ^5* in configuration 3d^, but there is such 
term in configuration Sd^As. Still, this assignment was 
used in the first calculation of Fe il [l3| for the identifica- 
tion of this level and a positive g-factor was obtained. In 
the next calculation Q this error was noticed and nega- 
tive g-factor for this level was reported. 

One of the main arguments which allowed to assign 
this level to the configuration Sd^AsAp was based on the 
analysis of f;-factors. The experimental g-factor for this 
level is 1.68 33], while all levels of the configuration 
3d^4p have significantly smaller (/-factors. On the other 
hand, the Lande value of the (7-factor for the level ^^7/2 
for the pure LS'-coupling scheme is 1.714. Calculated 
value of the g-factor appeared to be very close to the 
Lande value (see Table This was used in Ref. @ to 
estimate configuration mixing for this level to be about 
few percent and calculate corresponding correction to the 
q-factor: 



Sqg 



-180. 



(3) 



Since the publication of the paper Q the new NIST ta- 
bles for Fe II became available [33[ . These tables confirm 
assignment of the level ^Py2 62172 cm~^ to configura- 
tion Sd^AsAp. Moreover, the mixing of the configurations 
3d^AsAp and 3d^Ap for this level at a few percent was also 
independently confirmed from the analysis of the branch- 
ing ratios [3^. Therefore, we conclude that application 
of the correction ([3|) is justified. The correction of the 
opposite sign has to be applied to the closest interacting 
level 

At present it is impossible to reproduce this configu- 
rational mixing in the ab initio calculations because it is 
very sensitive to the relative position of the levels. The 
non-diagonal matrix elements which cause this mixing 
are very small, on the order of 10 cm~^, and mixing takes 
place only when the splitting between interacting levels 
is about 100 cm~^ or less. However, the accuracy of the 
present theory is insufhcient to reproduce the energy in- 
tervals between the levels of different configurations with 
such accuracy (see Table |l|. Because of that we still have 
to introduce correction ([3]) semi-empirically. If more ac- 
curate experimental data on g-factors and the branching 
ratios become available they may help to refine this anal- 
ysis. 



IV. CONCLUSIONS 

Our final results for the g-factors of the astrophysi- 
cally important levels of Fe 11 ion are given in Table IIVI 
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In addition to the large-scale 7-electron CI calculation 
for Dirac-Coulomb "no pair" Hamiltonian, we calculated 
corrections due to the Breit interaction and core-valence 
correlations. Both corrections appeared to be larger, 
than was anticipated in previous calculation [9]. Never- 
theless our calculations confirm previous conclusion that 
the g-factors of the levels of the configuration ScfiAp are 
large and positive, while the state ^Py2 62172 cni^^ 
belongs to the configuration SdPAsAp and has large neg- 
ative g-f actor. 

TABLE IV: Recommended values for g-factors in respect to 
the ground state, q^cc = gvai + Sqtr + Sqcv + 5qg, with un- 
certainties in parentheses (in cm~^). Valence contribution 
Qvai corresponds to the longest basis set [7spdf] from Table U] 
Breit correction Sqtr and core-valence correlation correction 
Sqcv are taken from Tables [IT] and IIIII Finally, g-factor cor- 
rection Sqg is given by Eq. ((3}. 

Qrcc 



Conf 


Level 


Qval 


Sqtr 


Sqcv 


Sqg 


This work 


Ref[9] 


3dMp 


^9/2 


1334 


19 


59 





1410 (60) 


1330 (150) 




6p)0 
^7/2 


1488 


11 


36 





1540 (40) 


1490 (150) 




^11/2 


1479 


11 


55 





1550 (60) 


1460 (150) 




^9/2 


1607 


-3 


55 





1660 (60) 


1590 (150) 


3d'^4p 


GjDO 
^7/2 


1134 


13 


394 





1540 (400) 


1210(150) 


3d^4p 


^7/2 


1240 


47 


449 


-180 


1560 (500) 


1100 (300) 


3d^4s4p 


6r>o 

^7/2 ~ 


-1506 


131 


168 


180 


-1030 (300) 


-1300 (300) 



We found out that the Breit interaction between va- 
lence electrons is of the same order of magnitude as the 
core-valence Breit interaction. Methods accounting for 
MBPT corrections developed so far do not allow us to 
treat core-valence correlations accurately for such com- 



plicated systems as Fe ii. Besides, we were unable to re- 
produce in the frame of the CI method strong mixing be- 
tween the level ^Py2 '^^ configuration 3d^4s4p and nearly 
located levels of configuration 3(i^4p. Corrections to q- 
factors due to this mixing were found semi-empirically 
using g- factor analysis This analysis is qualitatively 
confirmed by the analysis of the branching ratios [3J]. 
Considering the two problems mentioned above as a main 
source of uncertainties we can roughly estimate the final 
uncertainties. A, of the g-factors presented in Table HVl 
as 

A :^ ^{Sq^^r + 

As is seen from Table IIVI the uncertainties for the q- 
factors for the low-lying astrophysically important states 
arc at the level of 4% while for higher-lying levels they 
attain 30%. Further development of the CI-I-MBPT 
method allowing to treat core-valence correlations more 
accurately and appearance of more precise experimen- 
tal data on g-factors and transition rates for these levels 
would be very helpful in improving the accuracy of the 
q-factors. 
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